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1 Introduction 




is focused not on individual operators, but on paths in the space of such operators. We 
consider only paths connecting two operators conjugate by a scalar gauge transformation 
(so, they are loops up to a scalar gauge transformation). Such paths have a well known 
invariant, the spectral flow, which counts with signs the number of eigenvalues passing 
through zero from the start of the path to its end (the eigenvalues passing from negative 
values to positive one are counted with the plus sign, and egenvalues passing in the 
other direction are counted with the minus sign). The paper is devoted to the problem 
of computation of the spectral flow in the situation when all the operators along the path 
have the same symbol and the same boundary condition. 

Because these results are potentially useful for the physics of condensed matter, the 
author attempted to avoid advanced mathematical terminology and to explain the re- 
sults and the ideas behind their proofs in a way accesible to non-mathematicians. By 
the same reason, the author present the case of Dirac operators (Theorems [T] and |2) be- 
fore dealing with the more general case of Dirac type operators on domains equipped 
with an arbitrary metric. Note that physicists sometimes use more general boundary 
conditions than the ones considered in this paper. For example, the so-called armchair 
boundary conditions for graphen are of the type considered in this paper, but the zigzag 
boundary conditions for graphen are not. While it is not completely obvious, boundary 
conditions considered in this paper are just another forms of locally elliptic boundary 
conditions used in physics, as explained in Section HI Besides, we explain in Section 
|9] how the spectral flow can be computed in terms of general boundary problem for 
two-dimensional spinors considered by Akhmerov and Beenakker in HTJ. As illustration 
we show that the spectral flow vanishes if boundary condition does not break the time 
reversal symmetry. 

We start with the following situation. Let X be a compact planar domain bounded 
by m smooth curves (topologically it is a disk with m — 1 holes). Our operators act on 
spinor (i. e. spinor-valued) functions, which we identify with column vectors of two 
complex-valued functions: 



where q is a smooth function from X to C, x = (x^x 2 ) £ X, 3^ = 3/3x l . Our focus is on 
1-parameter families Dt of such operators parametrized by t £ [0,1]. In such a family 

1 In particular, boundary conditions defined by general pseudo-differential operators are not allowed. 




A Dirac operator acting on spinor functions has the form 
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the first term D involving derivatives is always the same, but in the second term the 
function q is allowed to continuously change with t, i. e. q = qt, where t E [0,1]. In 
agreement with the above, we assume that = |j.Do|J. 1 for some smooth scalar gauge 
transformation \i: X — > U(l). 

All operators Dt are considered with the same boundary condition of the form i(ui + 
vn.2)u + = B(x)u~, where n = (ni,U2) is the outward conormal to the boundary, and B 
is a real-valued smooth function on the boundary of X without zeros. Our first main 
result, Theorem [TJ asserts that the spectral flow of such a family of operators is equal to 
c m.YJ^=i bjl^j, where c m is an integer constant depending on m only, jXj is the degree of 
the restriction of \x to j-th connected boundary component, bj = 1 if B is negative on the 
j-th boundary component, and equal to otherwise. 

After considering this most special and very important situation, we turn our atten- 
tion to the situation of Dirac operators acting on N -dimensional spinor functions 



where, as before, X is a compact planar domain bounded by m smooth curves. A Dirac 
operator acting on N -dimensional spinor functions has the form 



In is N x N unit matrix, and QJx) is complex self-adjoint 2N x 2N matrix smoothly 
dependent on x e X. Again, our focus is on 1-parameter families D t of such operators 
parametrized by t E [0,1]. In such a family the first term D involving derivatives is 
always the same, but in the second term the matrix Q is allowed to continuously change 
with t, i. e. Q = Qt, where t e [0,1]. We assume that Di = ioDomT 1 for some smooth 
scalar gauge transformation \x: X — > U(l), where U(l) is considered as the subgroup of 
U(2N) consisting of the diagonal matrices with equal diagonal elements. 

All operators D t are considered with the same boundary condition i(ni + in2)u + = 
B (x)u , where B is a smooth map from the boundary to the space of complex self-adjoint 
invertible N x N matrices. Note that a local boundary condition is locally elliptic if and 
only if it can be written in such a form with B(x) invertible at any x; this boundary 
condition is self-adjoint if and only if B(x) is self-adjoint at any x. 

Our second main result, Theorem [2} asserts that the spectral flow of such a family of 
operators is equal to c m ^jl^ bj Uj, where c m is the same constant as in Theorem [1] (in 
particular, c m does not depend on the dimension N), [Xj is the degree of the restriction 
of \l to j-th boundary component (this restriction gives us the map from the circle to 
the circle because yi is a scalar gauge transformation), and bj is the number of negative 
eigenvalues of B(x) (counting with multiplicities) on the j-th boundary component. 




D 



D + Q(x), D 



i (cti3i + 0292) / 



where 
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Theorem [3] extends Theorem [2] to a still more general class of operators. These Dirac 
type operators involve in their definition an arbitrary (not necessarily flat) metric on 
X, and have the principal symbol defined by a Clifford multiplication which does not 
necessarily agree with this metric. While considering arbitrary metric is important for 
some physical applications, considering Clifford multiplication which does not agree 
with the metric on X does not seem neseccary. Nevertheless, we take care of this more 
general case because the proofs of our results crucially depend on its consideration. 
Moreover, we cannot prove Theorems [T] and |2] without proving Theorem [3] first. 

Notice that scalar gauge transformations \i: X — > U(l) leave invariant every local 
boundary condition, as well as the first term D of the operator D + Q(x): \x[D + QJ p. -1 = 
D + Q' for some function Q'(x). So any operator Do can be connected with the conjugate 
operator Di = |j.Do|J. 1 by the path (D + Qt(x)) with fixed boundary condition. On the 
contrary, non-scalar gauge transformations \x: X — > U(2N) do not have such properties, 
so the problem of computation of the spectral flow can not be stated in such a form as 
described above. If we allow general non-scalar gauge transformations, then we have to 
allow paths of operators (Dt) and of boundary conditions (Bt) with Bt and symbol of Dt 
being dependent on t. Some results about this more general case are outlined in Section 

m 

Note that in this paper the spectral flow is computed only up to multiplication by 
an integer constant c m depending only on m. For a disk with one hole (m = 2) the 
eigenvectors and hence the spectral flow are calculated explicitly in a special case; this 
is sufficient to determine C2; it turns out that C2 = 1 (see Theorem HJ. For the case m > 2 
this method fails because Fourier transform gives no help here. Nevertheless, the author 
expects that c m = 1 for all m; some reasons in favour of this conjecture are provided 
after proving that C2 = 1. 

Part I 

The spectral flow for Dirac operators 

2 The spectral flow 

Let H be a complex separable Hilbert space, (A t ), t e [0, 1] be a continuous 1-parameter 
family of bounded self-adjoint (or, what is the same, Hermitian) Fredholm operators in 
H. Near zero every A t has discrete real spectrum, which changes continuously with the 
variation of t. Hence one can count the net number of eigenvalues of A t passing through 
zero in positive direction as t runs from to 1, that is, the difference between the numbers 
of eigenvalues (counting multiplicities) crossing zero in positive and negative directions. 
This net number is called the spectral flow sf (A t ). The description of this notion can be 
found in EH. 

The case when Aq or Ai has zero eigenvalue requires some agreement on the count- 
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ing procedure; we use the following convention: take a small e > such that Ao, Ai have 
no eigenvalues in the interval [—£,0), and define the spectral flow as the net number of 
eigenvalues of A t + el which pass through zero. 

Let now (At) be an 1-parameter family of (not necessarily bounded) self-adjoint Fred- 
holm operators in H. For example, it can be a family of symmetric elliptic differential 
operators A t acting on sections of Hermitian bundle E over closed (that is, compact with- 
out boundary) manifold X. The definition of the spectral flow can be adjusted to this 
case, though more accurate consideration is needed, particularly due to the presence of 
various natural topologies on the space of such operators EUSHU. 

When a manifold has non-empty boundary, we have to consider the family (At, Bt), 
where At is a formally self-adjoint differential operator, and Bt is a "good" (self-adjoint 
elliptic) boundary condition for At at any t. One can see the notion of self-adjoint elliptic 
boundary value problem for operators of Dirac type in U3J. [71, and for general first order 
elliptic operators in [8j. 

Such differential operator At with boundary condition Bt defines the unbounded 
self-adjoint Fredholm operator on L 2 (X, E), which has unbounded discrete real spectrum. 
Intuitively, the spectrum of (A t , B t ) changes continuously with the variation of t, so the 
definition of the spectral flow works in this case as well EHH. However, the proof that 
the definition and the standard properties of the spectral flow are correct is considerably 
more difficult in this case. The crucial ingredient is the continuity (in t) of the family 
(A t , Bt) in the space of unbounded self-adjoint Fredholm operators on L 2 (X, E) with an 
appropriate metric. This was proved in (see [8], Theorem 7.16). This continuity 
property allows one to use the theory developed in HH5]] in full force. Our proof of 
Theorem |3] (see Part HT)) crucially depends on this theory, and, in particular, on Theorem 
7.16 from (8]|. The results of this theory needed for the proof of Theorem [3] are isolated 
in Section [TT] as properties (P0-P4). 

Note that if the spectra of (Ao, Bo) and (Ai,Bi) are the same (isospectral operators), 
which is the case in this paper, then there is another way to define the spectral flow of 
(A t , B t ). The set 

{(t,A) : A is the eigenvalue of (A t , B t )} 

can be uniquely represented as the union of the graphs of functions \{t) such that 
At(t) ^ Aj(t) for i ^ j. These functions give us the bijection (one-to-one correspondence) 
of the spectrum of (Ao, Bo) to the spectrum of (Ai, Bi). If these spectra coincide as 
subsets of R then this correspondence gives us the shift of the spectrum on the integer 
number of positions. This number is the spectral flow of (At, B t ). It is worth to note 
that for the isospectral case one can replace the level A = by any real number, and the 
difference between eigenvalues crossing the level in positive and negative directions will 
be the same |2]. 
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3 Dirac operators: the simplest case 

Suppose X is a compact planar domain bounded by m smooth curves (topologically it is 
a disk with m — 1 holes). We will use the notations x = (x 1 ,* 2 ) G X, 3^ = 3/3x\ 
Let us consider the Dirac operator on X 



(1) D = — i 



3i — t3 2 

ai + ia 2 o 



acting on a spinor function u: X — > C 2 , u 



A Dirac operator with non-zero vector potential has the form 

D=D + Q(x), where QM = ( q ° x) ^ 

q is a smooth function from X to C. 

Let X — y U(l) be a gauge transformation; we suppose that |x(x) G C, ||x(x)| = 1 for 
x G X. Let us take a Dirac operator Do = D + Qo(x) and connect it with the conjugate 
operator 

-l in , n , f i M- _1 (Sim- — i9 2 M-) 



by an one-parameter family of Dirac operators 



(2) D t = D + Q t , where Q t (x) 



q t (x; 
qtW 



qt is a smooth function from X to C continuously depending on t, t G [0, 1], 

qi-qo = i-M- -1 {d 1 [i + id 2 [i) . 
A self-adjoint locally ellipticJl local boundary condition for Dt has the form 

(3) in(x)u + = B(x)u~ on 3X, 

where B : 3X — > R \ {0} is a smooth function defining the boundary condition, n = 
(rii/Ti^) is the outward conormal to the boundary 3X of X at point x, and we identify n 
with the complex number + in 2 in ©. 

Note that tl\, n 2 coincide with the components of the outward normal to 3X for the 
case of Euclidean metric considered both here and in the next section. In Section [5] we 
consider a more general case of arbitrary metric on X, and the distinction between a 
normal and a conormal becomes essential there. 



2 Another name for "locally elliptic boundary condition" is "Sapiro-Lopatinskii boundary condition" 
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Remark. D + Q(x) is the Dirac operator on the trivial 2-dimensional complex vector 
bundle over X with compatible unitary connection defined by the function q(x). So the 
change of qt with t is equivalent to the change of the connection. 

Boundary condition © is gauge invariant with respect to the conjugation by \x, while 
Do and Di are conjugate by \i. So the operators Do, Di with the same boundary con- 
dition © are isospectral, and the spectral flow of the family D t gives us the shift of the 
spectrum of Dt when t runs from to 1. 

Theorem 1. The spectral flow of the family (D t ), t e [0, 1], with boundary condition © is equal 
to 

m 

3=1 

where c m is the integer constant depending on m only, [ij is the degree of the restriction of \x to 

f 1, ifB <0on 3Xj 
] ~ \ 0, ifB >0on 3Xj 

Here 3Xj are the connected components of the boundary of X, equipped with the orientation in 
such a way that the pair (outward normal to 3Xj, positive tangent vector to 3Xj) has the positive 
orientation on the plane (x^x 2 ). 




Figure 1: The case of two holes 

Note that since B ^ 0, it has definite sign at each boundary component 3Xj, so the 
constants bj are correctly defined. The restriction of \i to j-th connected component of 
3X gives us the map from the circle 3Xj to the circle U(l); |Xj is the degree of this map. 

This theorem follows from more general result which we formulate below. The gener- 
alization goes in two directions: (1) we allow arbitrary dimension of unknown complex 
functions u , u + , (2) we replace Dirac operator by operators of more general form. The 
value of C2 is computed in Section [6l 

Remark. Boundary condition © coincides with the boundary condition of Berry and 
Mondragon for the "neutrino billiard" [9j up to replacement of B by B _1 . In physical 
terms, one-parameter family of Dirac operators © describes the situation of continu- 
ously varying magnetic field so that the following two conditions are fullfiled: 
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(1) magnetic field at t = 1 coincides with magnetic field at t = all over the interior 
of X, 

(2) the fluxes through j-th hole at t = 1 and at t = differ by integer number (Xj in the 
units of the flux quantum. 

Let j = m corresponds to the outer boundary component and j = 1, . . . , m — 1 enumerate 

the holes. Considering that \x m = —YJ\=\^ M-j' we can reformulate Theorem [T] as follow- 
ing: the spectral flow of the operators family (O with boundary condition ® is equal 
to 

m-l 

Cm Y (bj-b m )Hj. 

5=1 

Thus the variation of magnetic field through j-th hole contributes to the value of the 
spectral flow with coefficient c m (bj — b m ). 

If the signs of B are the same on all boundary components, then the spectral flow is 
zero, no matter how magnetic field is varied (if only conditions (1-2) above are fullfiled). 
In the contrary, if B takes positive values on some boundary component and negative 
values on another, then we can vary magnetic field so that the spectral flow does not 
vanish. 



4 2N -dimensional Dirac operators 

Let X be as in the previous section. The standard 2N -dimensional Dirac operator has the 
form 

(4) D = -i [oidi + o 2 d 2 ) , where Oi = ^ *™ J , g 2 = 

In is N x N unit matrix. 

We will consider operators of the form D = D + Q(x) acting on spinor functions 

(5) u=^,u ± :X^C N , 

where Q is a smooth map from X to the space H(C 2N ) of complex self-adjoint (or, what 
is the same, Hermitian) 2N x 2N matrices. 

A self-adjoint local elliptic boundary condition for the operator O + Q has the form 

(6) in(x)u + = B(x)u~ on 3X, 

where B is a smooth map from 3X to the space of complex self-adjoint invertible N x N 
matrices, n = (111,112) is the outward conormal to 3X at point x, and we identify n with 
the complex number ni + rn.2. 
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The equivalent way of imposing boundary condition ^ is 



(7) (i(n 1 a 1 +n 2 a 2 )+^ B () 1 _° B j j u = an 3X. 

Theorem 2. Let Qt(x) be a continuous 1-parameter family of self-adjoint 2N x 2N matrices 
smoothly dependent on x e X such that D + Qx = \i (D + Qo) M- -1 for some smooth gauge 
transformation \x: X — > U(l). Let B be a smooth map from 3X to the space of complex self- 
adjoint invertible N x N matrices. Then the spectral flow of the family (D + Q t ) with boundary 
condition ((6]) is equal to 

m 

j=l 

ztfere c m zs the integer constant depending on m only, \Xj is the degree of the restriction of \x to 
the j-th connected component 3Xj of the boundary, bj is the number of negative eigenvalues o/B 
(counting multiplicities) on 3Xj (this number is correctly defined due to nondegeneracy o/B). 

This result is the corollary of Theorem |3] from the following section. 



5 Dirac type operators 

Let X be a compact planar domain bounded by m smooth curves and equipped with a 
Riemannian metric g (which is not necessarily flat). 

We call a first order formally self-adjoint operator D over X a Dirac type operator if its 
symbol has the form 

(8) 9 ={<>A=w<^ 



92/ V°2 

where (I) is a smooth map from X to the group GL + (2, R) of real invertible 2x2 matrices 
with positive determinant, and the matrices o\, 02 are defined by formula <(4]). 

In other words, Dirac type operator is the operator acting on spinor functions © and 
having the following form: 

(9) D = D 0/ Q = -i(pi(x)ai + p2(%)a 2 )+iR<DM + Q(x), 

where Q is a smooth map from X to the space H(C 2N ) of complex self-adjoint 2N x 2N 
matrices, 

1 r 1 

R<D = ^ [(Pl 9 l + P2 9 2) + (Pl3l + P23 2 ) 

(superscript t denotes the operation of taking the formal adjoint operator). More ex- 
plicitely, 

Ro(x) = ~ [9 a (Vg Pl ) + 3 2 (Vgp 2 )] g H(C 2N ), 
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where y/g = yj det(g^), the matrix (g^) is inverse to the matrix (g 1 ') = ((dx 1 , dx^) ), 
y/g dx : dx 2 is the volume element on X (of course, g^, g x \ and y/g depend on x). 

By D = T>x, g , n we denote the space of all operators having the form ((D for fixed X, 
g, N. Note that Dirac type operator (that is an element of ©x,g,N) is uniquely defined 
by the pair (<B,Q). 

A self-adjoint elliptic local boundary condition for Dirac type operator © has the 
form 

(10) in'(x)u + = B(x)u~ on 9X, 

where B is a smooth map from 3X to the space of complex self-adjoint invertible NxN 
matrices, the complex-valued function n' on 3X is defined by the formula n' = n[+ in 2 
with (nj,TL 2 ) = (111,112) O and n = (ni,U2) being the outward conormal to 3X at x G 3X. 
Recall that = Y_ 9^ f° r the components (rO) of the normal to the boundary 

Remark. Equation (fTD]> is just another form of the equation 

(11) ip + (x,n(x))u + = B(x)u~ on 3X, 
where 

PiX '^-^p+(x,£,) ) 

denotes the symbol £,ipi(x) + £,2P2M of the operator D in the direction of a covector 
£, = (£,1, £,2). Considering that in our case the operator p + (x, £,) is scalar, and p + (x, n(x)) = 
n'tx)!^, boundary condition (TTTJ may be written in simplified form (fTO)) . 

By ¥> = 2>x,n we denote the space of all smooth maps from 3X to the space of 
complex self-adjoint invertible NxN matrices. 

Suppose D e D, B 6 !B. We will write (D, B) for operator d9) acting on the domain 

ju G C 1 (X, C 2N ) : restriction of u to 3X satisfies boundary condition (fTO)) j , 

where C^X,^) is the space of continuously differentiable functions from X to C . 

Such operators have the following properties: 

1. For any D G T), B G 2 the operator (D, B) is (unbounded) essentially self-adjoint 
Fredholm operator, which has the discrete real spectrum. All its eigenvectors are 
smooth functions. (Lemma [H Section [TO) 
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2. Suppose QtM is continuous on (t,x), D G 2), B G 3. Then all the operators from 
the family (D + Q t , B) have the same domain, and this family is norm continuous 
in L 2 (X, g; <C 2N ). Therefore the spectral flow of the operators family (D + Q t , B) is 
well defined ([5], Proposition 2.2). 

Now we can formulate the main result of the present paper: 

Theorem 3. Let D G D be a Dirac type operator (|9]), B G 23 define boundary condition ((10) 
for D, QtM be a continuous 1-parameter family of self -adjoint 2N x 2N matrices smoothly 
dependent on x G X such that D + Qi = \i (D + Qq) yT 1 for some smooth gauge transformation 
[i: X U(l). TTien 

m 

Sf (D + Q t/ B) te[0;1] =c m ^bj^, 

j=i 

ztfere c m is the integer constant depending on m only, bj is ffoe number of negative eigenvalues 
o/B (counting multiplicities) on 3X,, ^ is the degree of the restriction of \x to 3Xj, 3X is oriented 
as described in the statement of Theorem [0 

Note that constant c m in all the Theorems [T]{3] is the same. 

Remark. Let S be a spinor bundle over X, ( ■ , ■ ) be an Hermitean metric on S compatible 
with its spinor structure, V be a connection on S compatible with its spinor structure 
and the Levi-Civita connection on TX. The Dirac operator on S in local coordinates has 
the form D = c(v) V v + c(w)V w , where (v, w) is a positive oriented orthonormal basis in 
T X X, and by c(v) we denote the action of a tangent vector v on spinors. 

The unitary skew-adjoint isomorphism J x = c(v)c(w) of S x does not depend on the 
choice of a basis (v, w) in T X X and defines the bundle decomposition S = S + © S , where 
S ± are the subbundles of S such that S x are the eigenspaces of J x corresponding to its 
eigenvalues =Fi. Due to the triviality of TX and of any complex bundle over X, we can 
fix some global positive oriented orthonormal basis field (v(x),w(x)) in TX and some 
trivialization of S . Let us extend the trivialization from S to S so that the action of the 
tangent vectors on the spinors in this trivialization has the form 

Then sections u of the spinor bundle S can be identified with the column vectors (0 of 
two functions : X — y C N , and Dirac operator D acting on such column vectors can be 
written in the form D = — i (piVi + P2V2), where pi, P2 are defined by formula ©, ®(x) 
is the transition matrix: (v,w) = (ei, e2)®(x), and by e\ we denote the vector (not the 
differential operator) 3^ to avoid misunderstanding. 

So any Dirac operator over X has the form ^\ with 0[x) satisfying the condition 
®(x)®*(x) = (g^(x)) and with the matrix Q(x) of very special kind. While considering 
arbitrary metric g and arbitrary Q(x) is important for some physical applications, con- 
sidering Clifford multiplication which does not agree with the metric on X (that is matrix 
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function <D(x) which does not satisfy the condition <D(x)<D*(x) = (g^(x))) does not seem 
neseccary. Nevertheless, we take care of this more general case because the proofs of our 
results crucially depend on its consideration. 



6 The case of one hole 

Here we compute the spectral flow for the case when X has just one hole (m = 2), and 
as a result find C2. 

Theorem 4. Q-i = 1. 

Proof. By Theorem |H the spectral flow does not depend on the geometry of X and on 
the choice of D e D, so it is sufficient to consider only the case when the computation 
is as simple as possible. Let us take the annulus X = {(r, cp) : 1 ^ r ^ 2} in the polar 
coordinates (r, cp) on the plane, with the metric ds 2 = dr 2 + d(p 2 , N = 1, 

n . / 3 r -i9<p\ up n ( it\ / +latr = l 

D = - 1 U + ^ J' ^ = e ' Q t= (^-it oJ' B = ( -latr = 2 

We obtain the following system for the eigenvector u and the eigenvalue A of (D + Qt, B): 

(— i9 r + 9(p — it)u+ = Au~ 
(— i9 r — 9cp + it)u~ = Au + 
u+ = iu _ at r = 1, 2 

All the eigenvectors of (D + Qt, B) are smooth functions, so we can seek them in the 
form u^r, cp) = Y-keZ, ^{r)e lk(p . Substituting it in the last system, we obtain 

9 r u^ — (k — t)u^ — iAu^ = 
9 r u k + (k — t)u k — iAu k = 
u k = iu k at r = 1, 2 

Equivalently, 

9 r (u+ + iu k ) = (k — t — A) (ut — iu k ) 
9 r (uj - iu k ) = (k - 1 + A) (u+ + iu k ) 
u k — iu k = at r = 1, 2 

and 9 2 (u k — iu k ) = ((k — t) 2 — A 2 ) (u k — iu k ) . So we have the following cases: 

• either u k = u k = 0, 

• or k — t + A = 0, u k = const, u k = iu k , 

• or (k - 1) 2 - A 2 = -(ttI) 2 , I e Z \ {0}, u+ - iu k = const • (e niXr - e~ nilr ^. 
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Therefore the set 



A = {(t, A) : A is the eigenvalue of (D + Q t , B)} 

is the union of the set Aj = {(t, A) : A — t e Z} (with the multiplicities 1 of the eigenvalues) 
and of the set A2 lying beyond the band |A| ^ n. 

If Aj (t) are the continuous functions from the interval [0, 1] to R such that A fl 
{0 ^ t ^ 1} is the union of the graphs of functions Aj (t) and At(t) ^ Aj (t) for i ^ j, 
then Aj (t) = j + 1 when —3 ^ j ^ 2 (up to shift of the numeration). So 

sf(D + Q t ,B) te[0/1] =1. 
On the other hand, by Theorem [3j 

sf (D + Q t/ B) te[0/1] = c 2 (bim + b 2 M-2) = c 2 [i 2 = c 2 

where by 3Xi, 8X2 we denote the inner and the outer boundary circles respectively 
Therefore C2 = 1. 

7 The case of several holes 

In this section we provide some evidence supporting the conjecture that c m = 1 for all 
m. 

Namely, let us realize X = X h as (m — 1) identical annuli arranged along the line and 
connected by the band of the width h, with the corners smoothed out, as on the Fig. |2l 




Figure 2: Contracting of the connecting band 



Let us consider the process of continuous decreasing of the band's width from h = 1 
to h = 0; we suppose that the annuli do not change in progress. Let us fix some function 
[i from X 1 to U(l) and take qt = it |a 1 (3i|x + i92M-)- Restricting \x and qt on X h , < h ^ 
1, we obtain operator |2]) over X h . Let us define the boundary condition by B h = +1 at 
the inner part U j<m 3X^ of 3X h and B h = -1 at the outer part 3X^ of 3X h . 

By Theorem [0 sf (D + Q t ) = c m ^j<m Mj does not depend on h. It is natural to 
suggest that the limit at h — > +0 of the (constant) spectral flow of the family (D + Q t , B h ) 
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for X h is equal to the spectral flow of (O + Qt, B°) for the "limit" domain X°, which is 
the disjoint union of m — 1 annuli, and the "limit" boundary condition B° = +1 at the 
inner boundary and B° = — 1 at the outer boundary of every annulus. 

However, sf (D + Qt, B°) for such union is the sum of sf (D + Q t , B°) for the annuli, 
and hence is equal to H, <m C2[ij = ci 2Ij< m V^- 

Therefore, if the assumption on the limit behavior of the spectral flow is true, then 

Cm = C2 a * an y m > 2- 



X 1 




Figure 3: Increasing of the holes 



Another way to have a look at the general case is to fix the outer boundary and to 
increase the holes up to their merging, as on Fig. |3] Here we obtain the single annulus in 
the limit of h = 0, and the same result c m = 1 if the passage to the limit will be justified. 

Alternatively, we can combine these two methods to obtain arbitrary number m/, 
1 ^ m/ ^ m — 1 of annuli in the end of the limit process, with the same result for c m . 

8 General case: first order elliptic operators 

The results of the present paper are concerned only with the case when X is a disk with 
holes. Easy modification of the proof gives us analogue of this result for the case of 
smooth compact oriented surface X with nonempty boundary, with the only change of 
c m to c m/ g, where c m/ g is the integer constant depending on the number m of boundary 
components of X and on the genus g of X. However, this still remains within the very 
restricted framework: all the operators D t are of Dirac type, both symbol of D t and 
boundary condition do not depend on t, conjugating gauge transformation is scalar. 

In fact, this result can be extended to much more general case. Namely, let X be a 
smooth compact surface, (At) be an 1-parameter family of first order symmetric elliptic 
differential operators acting on sections of unitary vector bundle E over X, and subbun- 
dle L t of E| ax defines a self-adjoint elliptic local boundary condition for A t at any t. 
Suppose that (Ai,Li) is conjugate to (Ao, Lo) by some gauge transformation \i (that is \x 
is unitary isomorphism of E, not necessarily scalar). Then operators (Aj, \-\), (Aq, Lq) are 
isospectral, and there arises the natural question about the spectral flow of the family 
(A t , L t ). This question will be considered in a forthcoming paper by the author IfTOll . In 
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that paper we will prove that 

m 

Sf (A t/ L t ) te [o,i] =c m , g Y_ <Pj/ 

j=l 

where c m , g is the integer depending on the number m of boundary components of X and 
on the genus g of X, cpj is the integer determined in a canonical way by the restrictions 
on j-th boundary component of the following data: 

(1) family (p t ), where p t is the symbol of A t ; 

(2) family (L t ) of boundary conditions; 

(3) gauge transformation p.. 

In particular, the spectral flow of (At, Lt) does not depend on the choice of the operators 
in the interior of X but only on the symbol of the operators on the boundary 

Theorem |3] of present paper fits into this general result as follows: c m = c m> Q, (pj = 
bj|j.j. Recall that p.j is invariant of the restriction of p to j-th boundary component of X, 
b, is defined from the restrictions of boundary condition and of the operator's symbol 
on j-th boundary component. 

9 The spectral flow for N = 2 in terms of the condensed 
matter physics 

In this section we compare our boundary condition ([6]) with "the general boundary 
conditions for the Dirac equation" for 4-dimensional Dirac operators formulated by 
Akhmerov and Beenakker in [lj. After that, we give some computations for the spectral 
flow in terms of (H. In particular, we show that the spectral flow vanishes in the case 
of time reversal symmetry (under the assumption of local ellipticity of the boundary 
problem). 

In this section we will temporarily use the notations from |1] in their original form 
and will formulate our results in the same terms. 

The long-wavelength and low-energy electronic excitations in graphene considered 
in (U are described by the Dirac equation HW = with Hamiltonian 

(12) H = vt ® [tr ■ p) 

acting on a four-component spinor wave function W = C^A/^b) (in our notations, W is 
two-dimensional spinor function, ^a = u+ / ^B = U- - , N = 2). Here v is the Fermi 
velocity, p = — ih-V is the momentum operator, cr ■ p = — ih (01V1 + 02V2), matrices 
%i, <Ji are Pauli matrices in valley space and sublattice space, respectively: 

0-0 = (J f),°l= (J J) ' 02 = (i 1 ) ' °" 3 = (0 -l) ' Ti = ai " 
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The general energy-independent boundary condition posed in [1] has the form 



(13) ¥ = on the boundary 

where M is a self-adjoint unitary 4x4 matrix depending on the point x G 9X and 
anticommuting with the current operator vto ® [cr ■ u-q). Here n-g is the outward normal 
to the boundary of X at x , so tib = (ui,U2) in our previous notations, and cr • n-g = 

+U202- 

Let us compare (|T3"1) with our boundary condition 

At first note that the condition "M is a self-adjoint unitary matrix anticommuting 
with the current operator" mean nothing but the condition of self-adjointness of the 
boundary problem (jl3|) . The authors of HI do not require local ellipticity from the 
boundary condition; however, in the absence of local ellipticity the spectrum of the 
operator is not expected to be discrete. Boundary condition (fl3]> is both locally elliptic 
and self-adjoint if and only if the matrix function M(x) can be represented by the formula 

for some complex self-adjoint invertible N x N matrix function B(x). For such M bound- 
ary condition (|T3~1) is equivalent to our boundary condition (O. 

The set of all possible self-adjoint unitary 4x4 matrices anticommuting with the 
current operator is paramertrised in [lj by the following 4-parameter family: 

(14) M = sin A t <g> [n\ ■ cr) + cos A {v ■ r) <g> [n 2 • cr), 

where A(x) is a "mixing angle", u[x), roi(x), n 2 (x) are unit vectors in R 3 = {(x^x 2 , x 3 )} 
such that n\ and n 2 are mutually orthogonal and also orthogonal to the boundary nor- 
mal tib(x), [u ■ t) = Y.i=i ylrv i> an d (^j • tr) are defined analogously Now we give the 



description of the ellipticity of the boundary problem (|13|) in terms of A, v, n\, n 2 , and 
compute bj as the functions of these parameters. 

From now on we will suppose that the frame [n%, n\, 112) is positive oriented in R 3 , that 
is its orientation coincides with the orientation of the frame [e\, Zi, e^) of basis coordinate 
vectors. This is possible because parameters {/K,n\,nj_,v) and (— A, — n\, n2, v) give us 
the same matrix M, so in a case of a negative oriented frame (n^, n\, n 2 ) we can change 
the signs of n\ and A simultaneously to obtain the positive orientation of the frame. 

Let cp(x) be a function from the boundary to the circle Rmod27i such that n 2 = 
sin cp • n + cos (p • e^, where e$ is the unit vector in R 3 in the direction of x 3 , r\{x) is the 
unit tangent vector to the boundary at x G 3X such that the pair (tib(x),t|(x)) has the 
positive orientation on the plane (x^x 2 ). Then n\ = cos (p • r\ — sincp • e$, and M is 
determined by the triple (A, cp, u). 

Proposition 1. Boundary condition (fi~3)) is locally elliptic for Dirac operator (|12)) if and only if 
A + (p ^ (mod 7t) and A — cp ^ (mod n) for any x G 3X. 



16 



In other words, the boundary condition is not locally elliptic if and only if ri 2 
± sin Ar] ± cos A e$ for some x G 3X and for some combination of signs ±. 



Proposition 2. If boundary condition (fl3|) is locally elliptic for Dirac operator (|12)) then it is 
equivalent to the boundary condition 

in(x)W A = B(x)¥ B 

with the matrix B(x) defined as following: 

„ „ „ „ , ^o±(£ / -t} „ A+cp „ A— cp 

B = (3+P+ + |3_P_, where P± = — — | (3+ = cot ^ , |3_ = tan — 

Here |3± are the eigenvalues ofB, P± are f/ze ortogonal projections on the invariant subspaces of 
B corresponding to the eigenvalues |3±. 

We prove these Propositions in the end of the section. 

Corollary 1. Let Qt(x) be a continuous 1-parameter family of self-adjoint 4x4 matrices 
smoothly dependent on x G X such that H + Qi = |j. (H + Qo) n _1 /or some smooth gauge 
transformation |x: X — >■ U(l). Suppose that boundary condition (fi~3)) is locally elliptic for Dirac 
operator ((T2|) . Then fhe spectral flow of the family (H + Q t ) wz'f/i f/zz's boundary condition is 
described by the formulae 



m 



Sf (H + Q t/ M) te[0 ,i] = c m Y_ b j h - ' 

i=i 

where c m , M-j are as in Theorem |2j bj depends only on the values of A, <p on j-th boundary 
component: 

{0, if both A + (p, A — (p belong to the interval (0, 7t) 
2, zf frof/z A + cp, A — cp belong to the interval {n, 2n) 
1, zjone o/ A + cp, A — cp belongs to the interval (0,7t) and another to the interval [n,2n) 

Proof. This follows immediately from Theorem |2] and Proposition |2] 

Let us inspect closer the case of time reversal symmetry The time reversal operator 
in the valley isotropic representation is 

T = -(T 2 ®ff 2 )e, 

with C the operator of complex conjugation [lj. The boundary condition preserves time 
reversal symmetry if M commutes with T. This implies that the mixing angle A = 
lU. By Proposition [1} in this case boundary problem (fi~3]> is locally elliptic if and only 
if n 2 (x) is not vertical for all x G 3X. If this is fullfilled then Corollary [T] allows us to 
compute the spectral flow regardless of other parameters: 
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Corollary 2 (the case of time reversal symmetry). Let Qt(x), t e [0,1] be a continuous 
1-parameter family of self -adjoint 2N x 2N matrices smoothly dependent on x e X suc/z £foz£ 
H + Qi = \x (H + Qo) [i^ 1 for some smooth gauge transformation \x: X — > U(l). Suppose that 
boundary condition is defined by formulas ((13)) , (fl4)> wzi/z A = 0, and that for any x 6 3X vector 
nzix) is not vertical. Then the spectral flow of the family (H + Q t ) is zero. 

Proof. By Corollary [TJ bj = 1 for all j. So we obtain 

m m 
Sf (H + Qt)te[0,l] =CmY_ b jh' = c m J" = 0- 

)=1 j=l 

Remark. Corollary |2] can be proved by other means as well, without use of Corollary [I] 
and of formula (|T4|) but using Theorem |2] directly. Namely, let M.W = W be locally elliptic 
boundary condition for Dirac operator ((12|) such that TMT 1 = M. At first note that 
the spectral flow of the family (H + Q t ) is independent of the choice of connection V, so 
we can assume that Vi = 3^. With this choice of connection, we have THT 1 = H. Let 
Q(( x ) =TQ t (x)T- 1 ,then 

Qi ~ Qo = T(Qi - Qo)^ 1 = TlyiHyT 1 - H)^ 1 = ^THT-V- THT 1 = ^'Hu/ _1 - H, 
where \i' = |a. _1 . By Theorem 12], 

sf(H + Q{,M) = c m ^b j( j j ' = -c m ^b j |i j =-sf(H + Q t ,M), 

where by sf (H + Q t , M) we denote the spectral flow of the family (H + Qt) with bound- 
ary condition NIW = W. If M commutes with T then 

(15) sf (T(H + Qt)T _1 ,TMT _1 ) = sf (H + Q{, M) = — sf (H + Q t , M). 

In Section [UJ we prove the conjugacy invariance of the spectral flow under unitary iso- 
morphisms of L 2 (X, g; C 2N ). Even though T is antilinear isomorphism of the Hilbert space 
L 2 (X, g;C 4 ), the spectral flow still remains invariant under conjugation by T. This can 
be proved using the uniqueness property of the spectral flow in the same manner as 
Property (P4) in Section [TTJ taking into account that conjugation by T preserves self- 
adjointness of operators in L 2 (X, g;C 4 ). Thus sf (T(H + Qt)T _1 ,TMT _1 ) coincides with 



Sf (H + Q t , M). Together with ©, this imply sf (H + Q t , M) = 0. 



Proof of Proposition \T\ By S = S(x) we denote the matrix [u ■ t) = Y.v=i ylx i- m our 
notations, 

, . . . / — sincp — iftcosaA . „ / cos cp — iftsincp\ 

M = sin A • t <g> . . M + cos A • S <g> . . I = 

\xn cos cp sm cp J \u\ sm cp — cos cp J 

f Si — inS2\ 
\inS2 —Si J ' 



18 



where Si = — sin cp sin A • I + cos cp cos A • S, S2 = cos cp sin A • I + sin cp cos A • S, I = I2 is 
2x2 identity matrix. Note that S 2 = I for any v, S 2 + S2 = I for any v, cp, A. 
Boundary condition NIW = W is equivalent to the following system: 



(16) 



-inS2U =(I — Si)u H 
inS2U+ = (I + Si)u" 



This boundary problem is locally elliptic for operator ((12)) if the linear space of the solu- 
tions of this system intersects with both spaces {u + = 0} and {u~ = 0} by zero subspace. 
This condition is equivalent to the invertibility of S2. Matrix S has the eigenvalues ±1, 
so S2 has the eigenvalues cos cp sin A ± sin cp cos A = sin(A ± cp). Both eigenvalues of S2 
are nonzero if and only if A ± cp ^ (mod n). This completes the proof. 

Proof of Proposition^ From (|T6|) we have B = S2 1 (I + Si). Taking into account the identity 
S 2 = I, we obtain 



S2 = (cos cpsin A — sin cpcos A) (cos cp sin A • I — sin cp cos A • S) , 

1 

S2 (I + Si) = ( cos cpsin A — sin cpcos A] (sin A — sin cp) (cos cp • I + cos A • S) 



sin A — sin cp 

cos cp • I + cos A • S) 



sin(A + cp) sin(A — cp 
Eigenvalues of S are ±1, so the eigenvalues of B are equal to 

sin A — sin cp 

(17) |3± = — — f -(cos cp ± cos A). 

sin(A+ cp) sin(A— cp) 

From last two formulas we have 

sin A — sin cp / , . , I + S , . . I — S 

B = . (cos cp + cos A) — h (cos cp — cos A) 



sin(A+ cp) sin(A- cp) \ 2 2 

= + = P+ P + + P- P - 

We can simplify (flT)) using sum-to-product trigonometric identities: 

sinA-sincp _ 2sin^ £ cos^ e ^_,„ A + ( P„_ A 

sin(A+cp)sin(A-cp) ~ (2 sin ^ cos ^) (2 sin ^ cos ^) ' 



2 sin — - — cos 



. n A+cp A-cp . A + cp . A-cp 

cos cp + cos A = 2 cos cos , cos cp — cos A = 2 sin sm . 

2 2 2 2 

Substituting this in ((17)) , we obtain 

A+cp „ A-cp 
(3 + = cot^^, (3_=tan^^. 
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This completes the proof. 



Part II 

Proof of Theorem SI 

Note that for D' = D + Q 0/ Q{ = Q t - Qo we have sf (D + Q t/ B) = sf (D' + Q(, B) with 
Qq = 0. By this reason, in the proof we will restrict ourselves by the families Qt with 
Qo = 0. 



10 Two technical lemmas 

First of all, we need to give some technical details. The reader interested only in the 
ideas behind the proof can go directly to the next section. 

Suppose D G D, B G 23. We will write (D, B) for the operator D acting on the domain 

(18) domain(D,B) = 

= ju G L 2 ^X;C 2N J : restriction of u to 3X satisfies boundary condition (flOl) j . 



Here L 2 (X;C 2N ) is the first Sobolev space; its elements are functions u G L 2 (X;C 2N ) 
such that 3iu, G L 2 (X;C 2N ). Strictly speaking, we use here not the restriction 
in the usual sense (trace map u i-> u| 9X ) but the extension by continuity of the trace 
map C°°(X;C 2N ) -> C°° (3X;C 2N ) to the bounded linear map from L 2 (X;C 2N ) to 
L 2 /2 (3X;C 2N ) ®. 

Note that the operator (D, B) defined here is the closure of the operator (D, B) defined 
in Section [5] (see B8]], Proposition 2.9). Using of non-closed operators in the first part of 
the paper is explained by our intention to avoid the introduction of Sobolev spaces and 
of the extension of the trace map as long as possible. Due to the following Lemma, these 
two definitions give us the operators with the same eigenvectors, so this slight abuse of 
notation does not cause any troubles. 

Lemma 1. For any D G T>, B G 23 the operator (D, B) is (unbounded) closed self-adjoint 
Fredholm operator on L 2 (X, g;C 2N ), which has the discrete real spectrum. Moreover, all its 
eigenvectors are smooth functions. 

Proof. Let B be a smooth function from 3X to GL(N,C). Then for any D G 2), A G C 
boundary condition (|T0|) satisfies the Sapiro-Lopatinskii condition for D — A: the inter- 
sections of the subspace {u: m'(x)u + = B(x)u~} c C 2N both with {u: u~ = 0} and with 
{u: u + = 0} are zero at any x G 3X. By Proposition 2.9 from (8), <[10]> is strongly regular 



20 



boundary condition for D, so all eigenvectors of (D, B) in L 2 (X, g; <C 2N ) are smooth func- 
tions. By the same Proposition, (D — A, B) is a closed Fredholm operator for any A G C, 
so the spectrum of (D,B) is discrete. 
For any u,w G L 2 (X, g;C 2N ) we have 



(Du,w) L 2 — (u,Dw} L 2 



((Du,w) - (u,Dw)) ^/g dx^dx 2 



IaJI 



'3i (v/gpi^w) + 3 2 ( v / 9P2U-,w)) dx dx' 



d ( (^/gpi^w) dx 2 - ( v /gp 2 u,w) dx 1 ) = -i y/g ( (piu,w) dx 2 - (p2U,w) dx 1 ) 
v J Jax v J 



ax 



((rtiP! +n 2 p2)u,w) ^ds 



ax 



) v^ds = 
((u _ ,vn/vv + ) — (vn/u + ,w~)) -y/gds, 



where ds is the length element on 3X. So for any u, w G domain(D, B) 



(Du,w) L 2 — (u,Dw) L 2 



ax 



(u" (B-B*) w -)Vgd S/ 



and the operator (D, B) is symmetric on L 2 (X, g;C 2N ) if and only if B(x) is self-adjoint 
at any x. 

Let now w G domain(D, B)*. By Proposition 2.9 from 0, domain(D, B)* is contained 
in L 2 (X, g; C 2N ), so we can use the computation above: 



(Du, w) L 2 — (u, Dw) L 2 



u 



ax 



Bw )) v^ds 



for any u G domain(D,B). Therefore, in / w + — Bw~| ax = for any w G domain(D, B)*, 
domain(D,B)* = domain(D,B), and (D,B) is self-adjoint on L 2 (X, g;C 2N ). All eigen- 
values of a self-adjoint operator are real. This completes the proof. 

In the statement of Theorem [3] we used only norm continuous paths of operators 
with fixed domain. But for the proof of Theorem [3] we have to deal with the paths of 
more general kind, when neither symbol of the operator nor boundary condition are 
fixed any more. The paths we need for the proof are not norm continuous but only 
graph continuous (note that by Proposition 2.2 from [5] any norm continuous path is 
graph continuous as well). So further we will use the graph topology on the space of 
closed densely defined self-adjoint operators on a separable Hilbert space H (in our case 
H = L 2 (X,g;C 2N )). 

There are various definitions of the graph distance, all of which give the same graph 
topology 0]. One can take d G {A,A') = ||(A + iI) _1 - (A ; + il)~ 1 ||, or alternatively 
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cIg(A, A') = ||P A — P A /||, where Pa, Pa' are the orthogonal projections of H x H onto 
the graphs of A, A' respectively. 

Let us introduce the following metrics in D and I>: 

dfDo^D^q,) = ||Q-Q / || C (x) + ll (I) - 0/ |lci(X) = 

= max ||Q(x) - Q'{x)\\ +max (II® (x) - <D'(x)|| + ||3i<D(x) - 9i<D'(x)|| + ||3 2 0(x) - 3 2 / (x)||) , 

xeX xeX 

d(B / B / ) = ||B-B'|| cl(ax] =max(||B(x)-B / (x)|| + ||3 s B(x)-3 s B / (x)||), 

where s is the length parameter on 3X. Here we use any of the standard norms on the 
spaces B(C N ) and B(C 2N ) of complex N x N and 2N x 2N matrices, and on the space 
B(R 2 ) of real 2x2 matrices. 

Note that (Do t Q t ,B t ) is the continuous path in D x 23 if and only if QtM, ®tM/ 
B t (x), and the first partial derivatives of ® t (x), B t (x) with respect to x are continuous 
functions of (t,x). 

Denote by HF(H) the space of closed self-adjoint (or, what is the same, Hermitian) 
Fredholm operators in separable Hilbert space H. We fix graph topology on HF(H). 
Nevertheless we will usually write "graph continuous" instead of mere "continuous" 
for the maps to HF(H) to avoid a misunderstanding. 

By Lemma CO we have the natural inclusion DxS <— > HF (L 2 (X, g;C 2N )), which 
carries a pair (D, B) G T> x 23 to the operator D with the domain (|T8|) . 

Lemma 2. The natural inclusion BxB^ HF (L 2 (X, g;C 2N )) is graph continuous. 

Therefore, if t i — >■ (Dt, Bt) is the continuous path in D x 23, then (Dt, Bt) defines 
the graph continuous path in HF (L 2 (X, g;C 2N )), and the spectral flow of the operators 
family (Dt, Bt) is well defined. 

Proof. Let us consider the smooth map 

*: B(C«) H(C-), A~P=(_fc A ". A A ) C N I^A-a)" 

which carries A G B(C N ) into the orthogonal projection P of C 2N with KerP = 
{u = (u + ,u~) : u + ,u~ G C N ,u + = Au~}. It induces the continuous map 

ip*: C 1 (3X,B(C N )) ->• C 1 (3X,H(C 2N )) . 

Composing ip* with the continuous map 

D x S -> C 1 (3X,B(C N )), (D,B) ^ -ip + (x,n(x)) _1 B(x), 

we obtain the continuous map 

W: V x 23 C 1 (3X,H(C 2N )) , 
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which carries (D, B) into the orthogonal projection P of L 2 (3X, g| 9X ; C 2N ) with the kernel 
defined by boundary condition ((10|) 0. 

By Proposition II. 1.1 from ffTTfl , we have the continuous inclusion of the Banach spaces 

B(L? (3X;C 2N ))cB(L? /2 (3X;C 2N )), 

where B(V) denotes the space of bounded linear operators on a Banach space V, L 2 is 
the (fractional) Sobolev space. Composing it with the natural continuous inclusion 

C 1 (3X,B(C 2N )) C B (L 2 (3X;C 2N )) , 

we obtain that the map ^:Dx3->B (l 2 /2 (3X;C 2N )) is continuous. 

The natural map from D to the space of bounded linear operators from L 2 (X;C 2N ) 
to L 2 (X;C 2N ) is continuous too: 

||D<d,q - Do',Q'|| 1;0 < const M|® - <D'|| C(X) + ||R<d - R<D'|lc(X) + || Q - Q'llc(X)) ^ 

< const (||<D - O'HcipQ + ||Q - Q'Hccx)) ■ 

By Theorem 7.16 from [8J and by Lemma [H this implies that the inclusion !Dx!B4 
HF (L 2 (X, g;C 2 )) is graph continuous. This completes the proof. 



11 Basic properties of the spectral flow 

There can be different versions of the definition of the spectral flow when one or both of 
the endpoints of the path is non-invertible. If a path is a loop up to a gauge transforma- 
tion as in the first part of the paper, then the value of the spectral flow is independent of 
the choice of the definition. But for the proofs below we have to fix some choice. 

Definition. Let (At) be an 1-parameter graph continuous family of closed self-adjoint 
Fredholm operators in a separable complex Hilbert space H. Take a small e > such 
that Ao, Ai have no eigenvalues in the interval [—£,0). We put sf (At) be equal to sf (At + 
el), were we use any of the (equivalent) definitions of the spectral flow for the path of 
operators with invertible endpoints from |H [5]]. This definition does not depend on the 
choice of such e. 

We will need the following properties of the spectral flow. 

(P0) Zero crossing. In the absence of zero crossing the spectral flow vanishes. More pre- 
cisely, suppose y: [0, 1] — >• D x 25 is the continuous path such that is not the eigenvalue 
of y(t) for any t e [0, 1]. Then sf (y) = 0. 

3 Here we use the general formula for the orthogonal projection P with the kernel {u + = Au~} for 
arbitrary matrix A. Actually, in our case A = (in') _1 B is normal: AA* = A* A. 
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(PI) Homotopy invariance. The spectral flow along the continuous path y : [0, 1] — >■ 
BxS does not change when y changes continuously in the space of paths in D x 23 
with the fixed endpoints (the same as the endpoints of y). 

In other words, for the continuous map h: [0,1] x [0,1] 4 D x 2 such that H s (0) = 
(D , B ), Ml) = (Di,B!),wehavesf (ho(t)) te[0/ i] = sf [^i[t)) tem . 

(P2) Path additivity. Suppose y : [a, c] — > D x 23 is a continuous path, a ^ b ^ c. Then 
sf(y(t)) 

te[a,c] — S f h / ( t ))te[a,b] + Sf (y(t)) te [b,c]- 

(P3) Additivity with respect to direct sum. Let Ni, N2 be natural numbers, (D{, B|) be 
continuous paths in T>^ i x 23jsj i . Then the spectral flow along the path (D| © D 2 , B| © B 2 ) 
is equal to the sum of the spectral flows along the paths (D\, B|) and (D 2 , B 2 ). 

(P4) Conjugacy invariance. Let J±: X — > U(N) be unitary N x N matrices smoothly 
dependent on x <E X, J = P+ ®j : X ->■ U(2N), (D t/ B t ) be a smooth path in D x 23. 

Then Sf (D t , B t ) = sf (JDtJ- 1 , J-BtJI 1 ). 

More generally, if H is a separable complex Hilbert space, J is an unitary isomorphism 
of H, (At) is an 1-parameter graph continuous family of closed self-adjoint Fredholm 
operators, then sf (A t ) = sf (JA t J _1 ). 

Remark. Properties (PI) and (P2) imply that the spectral flow along the path is opposite 
to the spectral flow along the same path passing in the opposite direction. 

Proof. By Lemmas HHH the inclusion of D x 13 into HF (L 2 (X, g;C 2N )) is graph contin- 
uous. So it is sufficient to prove Properties (P0-P4) for graph continuous paths in the 
space HF(H) of all closed self-adjoint Fredholm operators in separable Hilbert space H; 
this will imply properties (P0-P4) for the paths in D x 13. 

First three properties of the spectral flow for graph continuous paths in HF(H) are 
proved in |H (Proposition 2.2), taking into account the convention from Section |2] for the 
case when y(0) or y(l) are non-invertible. 

Conjugacy invariance of the spectral flow for graph continuous paths in HF(H) fol- 
lows from the uniqueness property of the spectral flow. Namely, let J be an unitary 
isomorphism of a separable complex Hilbert space H. To each graph continuous path 
(A t ) in HF(H) assign the integer sf new (At) = sf(JA t J _1 ). Then sf new satisfies Concate- 
nation, Homotopy and Normalization properties in the sense of 0. By Theorem 5.9 
from [5], this imply that sf new equals sf for the paths in HF(H) with invertible endpoints. 
Taking into account our convention from Section |2] and choosing a small £ > such that 
Ao, Ai have no eigenvalues in the interval [— e,0), we obtain 

Sf (JAtJ- 1 ) = Sf (JAtJ- 1 + el) = sfnew(A t + el) = Sf (A t + el) = sf (A t ). 

To prove (P3), consider graph continuous paths (At), (A{) in HF(H), HF(H') respec- 
tively. Suppose for a while that Ao, Ai, Aq, and A[ are invertible. The path (A t © A() te [o,i] 
is homotopic to the concatenation of paths (A t © Ao)te[o,l] an d (Aj © A() te [ 0/1 ] in HF(H © 
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H'). The spectral flow of the path (A t © Aq) in HF(H © H') considered as the function of 
(A t ) satisfies Concatenation, Homotopy and Normalization properties in the sense of [5], 
so by the uniqueness property of the spectral flow from [5] we have sf (A t © Aq) = sf (At). 
Similarly sf (Ai © A[) = sf (A(). Therefore, sf (A t © A{) = sf (A t ) + sf (A{) for any paths 
(At), (A{) with invertible endpoints. Taking into account our convention from Section 
we obtain that sf(A t ©A() = sf (A t ) + sf (A{) for arbitrary paths (A t ), (A{). This 
completes the proof. 

12 Independence of the choice of family (Q t ) 

Let us prove that the spectral flow along (D + Qt, B) does not depend on the choice of 
(Qt) when D, B, \i are fixed. 

Let Qt, Q( be two 1-parameter families of smooth maps from X to H(C 2N ) such that 
Qo = Q >0 / Q 1 = Q{ = uT^- 1 -D. 

The path D + Q t can be continuously changed to the path D + Q{ in the class of 
paths in D with the fixed endpoints. For example, we can take the homotopy h(s, t) = 
D + (1 — s)Q t + sQ{. By the homotopy in variance property (PI) of the spectral flow, 
sf (D + Q t , B) te [o,i] = sf (D + Q(, B) te [o,i]. 

Therefore, if Q = 0, Qi = \xDyr 1 - D then 

Sf (D + Q t , B) te[0 ,i] = F(X, g, N, D, B, n) 

for some integer-valued function F. Now we will investigate the properties of this func- 
tion. 

13 Independence of the choice of operator D 

1. Suppose that Do is homotopic to Di in D, that is there exist a continuous 1-parameter 
family of the Dirac type operators D s connecting Do with Dj. We will show now that 
F(X, g, N, Do, B, |x) = F(X, g, N, Di, B, |x). 

Let us consider 2-parameter family of the Dirac type operators D s/t = (1 — t)D s + 
t[xD s |x : . Note that D s , = D s , D s4 = nDs^r 1 , D S/t - D S/0 = tQ s , where Q s = mD^ 1 - 
D s is the 1-parameter family of 2N x 2N self-adjoint complex matrices smoothly depen- 
dent on x G X. 

Let us define the path y\ : [0, 3] — > D by the formula 

f D t/0 , tG [0,1] 
YlW = { D u -i, tG [1,2] 
{ D 3 _ t/1 , tG [2,3] 

In other words, we consequently go around the left, top and right sides of the rectangle 
on Fig. H] in clockwise direction. The path yi can be continuously deformed to the path 
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Figure 4: Homotopy from yo(t) to yi(t) 



Do,t = Dq + tQo within the rectangle. For example, we can take as such a deformation 
the family 



On the other hand, the spectral flows along the first and the third parts of yj are mutually 
reduced by (P4): 

sf (y(t), B) te[0/1] + sf (y(t), B) te[2/3] = sf (D S/ B) s€[0/1] - sf {[iD s [r\ B) s€[0/1] = 0. 

Therefore, sf (D + tQ ,B) t€[0/1] = sf (yi(t),B) te[12] = sf (D x + tQi, B) t€[0/1]/ and 
F(X, g, N, D , B, n) = F(X, g, N, D l7 B, ^i). 

2. Now we will simplify D step by step. 

At first, we can continuously change D = Do,q = — i (Pi9i + P2^2) + t-Rcd( x ) + Q(^) 
to the operator Dq q, for example, along the path D^ ^^q. 

Further, take a smooth map h: [0, 1] x GL + (2,R) ->■ GL+(2,R) such that h(0, •) is the 
identity map, while the image of H(l, •) is the group SO(2, R) of 2 x 2 orthogonal real 
matrices with determinant equal to one (the existence of such a family is well known in 
any dimension; this is an easy application of Gram-Schmidt orthonormalisation proce- 
dure). The operator D$ ; o can be continuously changed in D along the path D h ( t ^ /0 to 
the operator where 
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for some smooth function cp from X to S . So we have 

F(X, g, N, Do 7 q, B, \i) = F(X, g, N, D o ,, , B, yi). 

On the other hand, D<d/ can be represented as J^D^qJ, where I = I2 is the identity 
2x2 matrix, 

n^Co + ;>( el o lN O eu(2N) - 

Let Qt be an 1-parameter family of self-adjoint 2N x 2N complex matrices such that 
Qo = 0, Qi = (j-Dfl^oM- -1 — D®',0- Applying property (P4) of the spectral flow, we obtain 

F(X, g, N, D o ,, , B, \i) = sf (D o , + Qt, B) = sf (J (D o , /0 + Qt) J~\ J BJZ 1 ) = 

= Sf (Di /0 + JQtJ -1 , B) = F(X, g, N, D I/0 , B, yi), 

because JQiJ 1 = \i (jD o / /0 J _1 ) M-" 1 - JD^oJ- 1 = uL^oH^ 1 - D I/0 . 

Therefore, F(X, g, N, D, B, \x) does not depend on the choice of D £ D, so from now 
on we will write F(X, g, N, B, \i) instead of F(X, g, N, D, B, 

14 Independence of the metric and invariance under the 
change of variables 

We prove here that F(X, g, N,B, \i] is independent from the choice of the metric g on 
X, invariant under the change of variables, and does not depend on the geometry of 
X, using the fact that the number of holes is the only topological invariant of the disk 
with holes, and that the spectral flow is conjugacy invariant and does not depend on the 
choice of the operator. 

Let X, X' be compact planar domains, each bounded by m smooth curves, and g, g' 
be Riemannian metrics on X, X' respectively 

As well known, there exists an orientation-preserving diffeomorphism f : X' — > X. □ 

We define 9 as the smooth function from X' to R + such that f * dvol = 9 dvol', where 
dvol, dvol' are volume elements on X, X' respectively. 

Diffeomorphism f defines the unitary isomorphism J from the Hilbert space 
L 2 (X, g;C 2N ) to the Hilbert space L 2 (X ; , g';C 2N ), u H- v^fu. | 

Isomorphism J transforms the operator D e T>x, g , n with symbol p to the symmetric 
operator D' = JDJ -1 on X' with symbol p'. For any x' E X', x = f(x'), any cotangent 

4 In other words, there exist a smooth one-to-one change of variables (x' 1 ,*' 2 ) =x'- l >x= (x^x 2 ) with 
the smooth inverse and with positive Jacobian determinant det(3x/3x'), which transforms X' onto X. 
5 As usual, by f * we denote the homomorphism from the differential forms (in particular, functions) on 

X to the differential forms on X', which is induced by f. In coordinate form, 6(x'l = detf 

vV( x ') 

6 That is (Ju)(x') = vWIu(f(x')). 
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vector £, £ T X *X, £/ = f*£„ we have pV/£') = p[x,l) f that is p'(x') = p(x) = 

(w) m coor diriate representation. The matrix (j^j is contained in GL + (2,IR) 

for any x G X, so D' G £>x',g',N- 

Let p. be a smooth function from X to U(l). Taking the map \i' = f*p from X' to U(l) 

and the map B' = \\f*n\\~} f*B from 3X' to H(C N ), we obtain 

p'D V'" 1 — D' = p/ (jDJ- 1 ) p'" 1 - IDT 1 = J (pDp" 1 - D) J" 1 . 

50 if Q t connects Q = with Qi = pDp _1 - D, then Q( = JQtJ 1 connects Q' Q = with 
Qj = (j.'D'(j. — D ; , and by the conjugacy invariance of the spectral flow (P4), we have 

Sf (D + Q t/ B) = Sf (j(D + Qtjr 1 , B') = Sf (D ; + Q[, B ; ). 

However, B' is homotopic to f*B in ®x',N/ while the spectral flow of (D' + Q(, B) is 
invariant under the continuous change of B in 23x' 7 n (this is verified in a way similar 
to the proof in Section [HI). Therefore sf (D ; + Q{, B') = sf (D ; + Q{, f*B), and finally we 
obtain 

(19) F(X, g, N, B, p.) = F(X', g', N, f*B, f». 

This completes the proof. 

In particular, for any two metrics g, g' on the same X, using the identity diffeomor- 
phism f, we have 

F(X,g,N,B, M .) = F(X,g , ,N,B, M .). 
Further we will write F(N, B, \x) instead of F(X, g, N, B, p.). 

15 Boundary conditions 

Let us investigate the dependence of F(N, B, \i) on B. 

F(N,B, \i) does not change when B continuously changes in !B; this is verified in a 
way similar to the proof in Section [13J 

Let bj be the number of negative eigenvalues of B (counting multiplicities) on 3Xj. We 
prove that the ordered set b = (bj). 1 uniquely determines B up to continuous variation 
of B in S. 

Obviously, b is invariant with respect to such variations, so we only have to prove that 
any two B, B ; with the same b are homotopic. It is sufficient to prove that any smooth 
map A from the circle S 1 to the space of complex self-adjoint invertible N x N matrices 
is homotopic (in the space of all such maps with C 1 -metric) to the steady map sending 

5 1 to the point (—Ik) © iN-k £ H(C N ), where k is the number of negative eigenvalues of 
A(x), x G S 1 . 
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1. Let us consider the continuous 1-parameter family A s of smooth maps from S 1 to 
the space of complex self-adjoint invertible NxN matrices defined by the formula A s = 

A • ((1 — s)In + sA 2 ) . This expression is correct because (1 — s)Im + sA 2 is self-adjoint 
and positive definite for any s G [0, 1] . The family A s gives us the deformation from 
A = Aq to smooth map Ai from S 1 to the space of self-adjoint unitary NxN matrices. 

2. The connected component of Ai (x) in the space of self-adjoint unitary NxN matrices 
is diffeomorphic to the space Grc(k, N) of all k-dimensional linear subspaces of C N . This 
diffeomorphism is defined by the correspondence U h-> Ker(L\j +U), which associates 
with U the invariant subspace V C C N of U corresponding to eigenvalue — 1 of U. The 
inverse diffeomorphism is defined by the formula Vh-U = (— I)v © ly- 1 -- 

The complex Grassmanian Gr<;(k, N) is known to be simply connected, so any two 
continuous maps from the circle to Gr^ (k, N ) are homotopic. Taking into account that 
Grc(k, N) is the smooth manifold, we obtain that the space of smooth maps from the 
circle to Gr^ (k, N ) (with C 1 -metric) is path-connected. The same is true for the connected 
component of the space of self-adjoint unitary NxN matrices which is diffeomorphic to 
Gr<c(k, N), so A can be continuously changed in the class of smooth maps to the steady 
map x !->■ (—lie) © In-1c- This completes the proof. 

16 Gauge transformations 

1. We will prove that F is linear in \i, that is F(N, B, \i\\i2) = F(N, B, m) + F(N, B, \i 2 ) for 
any smooth functions \i2'. X — >■ U(l). 

Let Qi = (j-tD^r 1 — D. Then Qi + Q2 = (M-1M-2) D {\i\\x-2)~ — D, so by definition 
F(N, B, is equal to the spectral flow along the path (D + P t , B) te [ 02 ], where Po = 0, 
P2 = Qi + Q2- We can take P t composed from two parts: from to Qi and then from Qi 
to Qi + Q2, for example, 

p = / tQi, te[0,l] 
x \ Qi + (t-l)Q 2/ te [1,2] 

Using the property (P2) of the spectral flow, we obtain 

F(N, B, (im 2 ) = Sf (D + P t , B) te[0/1] + Sf (D + P t , B) te[1/2] = 

= Sf (D + tQx, B) t6[0/1] + Sf ( (D + Qi) + tQ 2 , B) t6[0/1] = 

= F(N,B,m) + F(N,B,^ 2 ), 

so F is linear in \i. 

2. By M denote the set of equivalence classes of smooth functions |x: X — > U(l), where 
two functions are equivalent if one of them can be continuously changed to another in 
the space of smooth functions from X to U(l) (with C 1 -metric). We will consider M as 
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the Abelian group, where the group structure on M is induced by the group structure 
onll(l). It is well known that 

M = {(n 1/ .. ; ^)eZ m : H^j = }/ 

with the group structure induced from Z m , and the class of (J. in M is defined by the 
m- tuple p. = (\Xj), where \ij is the degree of the restriction of (x to 3Xj. 

Let us prove that F(N, B, \i) depends only on the class of \i in M. 

Suppose that \±t is a continuous path in the space of smooth functions from X to U(l) 
such that (J-o(x) = 1. By the previous clause, it is sufficient to prove that F(N, B, = 0. 
Let us take Q t = [XtDM-t 1 — D. Taking into account that Qi = (j-iD^ 1 — D, we obtain 
F(N,B, = sf (D + Qt,B). But all the operators (D + Q t/ B) are conjugate to (D,B) by 
(j-t and therefore are isospectral. Let e > be such that (D,B) has no zero eigenvalues 
in the interval [-£, 0). Then sf (D + Q t , B) = sf (D + Q t + eI 2 N/ B) = by (P0) because all 
the operators (D + Q t + £i2N/ B) have no zero eigenvalues. This completes the proof. 



17 Bilinearity 

In the previous sections we have proven that F depends only on the integer numbers N, 
bi, . . . , b m , (J-i, . . . , |x m . Now we will study this dependence more closely. 
By S denote the set of all possible (m + l)-tuples (N, bi, . . . , b m ): 

S = |(N,b 1 ,...,b m ) e Z m+1 : N ^ l,0Oj ^ n}. 

F defines the map from S x M to Z (which we denote by the same letter F for simplicity) 
satisfying the following conditions: 

F (N, b, p.© p/) = F (N, b, p.) + F (N, b, fi') 
F(N + N / ,b©b / ,fl) =F(N,b,fl) +F(N / ,b / ,fl) 

where p. = . 1 m , b = (bj). 1 m , symbol © denotes the componentwise addition. 
Indeed, the first equality has been proven in section [T6]. while the second equality is by 
the property (P3) of the spectral flow. 

Hence F is a bilinear function, and therefore there is a homomorphism from Z m+1 © 
M to Z such that F can be represented as the composition 

(20) SxM4 Z m+1 x M — > Z m+1 © M Z, 

where the first arrow is induced by the natural embedding of S into Z m+1 , and the 
second arrow is the canonical map of the direct product to the tensor product. 
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Let us consider operator (O with boundary condition If (D + Q t ) u = and 
i (ni + XTI2) u + = Bu on 3X, then 



(B(x)u ,u ) ds 

9X 



ax 



'i (ui + rn.2) u + ,u ) ds = 



((— i (82 + 182) + qt) u + ,u )dx x dx 2 — (u + , (— i (3j — 182) + q t ) u )dx 1 dx 2 = 0, 

Jx 

where ds is the length element on 3X. 

Suppose now that the sign of B is the same on all boundary components. Then 
from the last equality we have vT = on 3X, u + = —\{n\ — vn.2)Bu~ = on 3X. 
Thus u = on X by the weak inner unique continuation property of Dirac operator 
0. So (D + Qt, B) has no zero eigenvalues at any t for such B, and by Property (PO) 
sf (D + Qt, B) = 0. Finally we obtain F(l, 0, p.) = F(l, 1, p.) = at any p., where we denote 

o = (o,...,o), i = (i,... / i)GZ m . 

Let us consider the group M/ which is quotient of Z m+1 by subgroup spanned by 
elements (1,0) , (l, l) £ Z m+1 . Note that M' coincides with the quotient group Z m / (l), 
so it is naturally isomorphic to the Abelian group Hom(M, Z) of all homomorphisms of 
MtoZ. 

By previous arguments, there exists such homomorphism F: M/ ® M — > Z that F is 
the composition of the following homomorphisms: 

(21) S x M •— > Z m+1 x M — > Z m+1 ® M M' g> M 4 Z, 

where the first two arrows are the same as in (|20l> , and the third arrow is induced by the 
natural projection Z m+1 — > Nl'. 



18 Invariance under the action of symmetric group 

Let Diff + (X) be the group of all diffeomorphisms of X preserving orientation, f £ 
Diff+(X). As it was shown in Section HI F (N, f*B, f» = F(N, B, and hence 

F(N,f*b\f*fl) =F(N,b,ft), 

where f * acts on b and p. by the permutation of the coordinates, corresponding to the 
permutation of the boundary components of X by f. It is well known that any permu- 
tation of the boundary components of X is realized by some element of Diff + (X). Thus 
F(N,b, p.) is invariant under the action of symmetric group S m (the group of permuta- 
tions of ra elements) on (b, ft) by the permutations of the coordinates. 

All permutations of the coordinates leave he element 1 of Z m invariant, so S m acts 
onM' = Z m / (l) in exactly the same way, and F is invariant under the action of S m , too. 

Extending F by linearity from M ; <g> M to V <8> V, V = M' ® C = C m /(l), V = 
M C = {v £ C m : = 0}, we obtain homomorphism F^ : V <8> V — > C, coinciding 
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with F on the lattice M/ <g> M c V ® V. Obviously, is invariant with respect to the 
action of S m on V ® V as well. 

V and Hom c (V,C) coincide as the representations of S m , so the vector space of all 
invariant homomorphisms from V <g> V to C is isomorphic to the vector space of all 
equivariant homomorphisms V — > V. But the latter space is 1-dimensional by Schur's 
lemma, because V is the irreducible representation of S m | fl2| . So F c (v ; <g> v) = c Y_) v \ v ) 
for some constant c G C, and F(N, b, p.) = c Y_ bj(^, where c depends only on m. 

On the other hand, F is integer-valued and, in particular, c = F (1, (0, 1), (—1,1)) e Z. 

Finally, we obtain sf (D + Q t , B) te [o / i] = c m YJ\Li ^°)V i )> where c m is the integer constant 
depending on m only, and Theorem [3] is proved. 
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